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Abstract — This paper proposes a control allocation technique 
that can help pilots recover from pilot induced oscillations 
(PIO). When actuators are rate-saturated due to aggressive pilot 
commands, high gain flight control systems or some anomaly 
in the system, the effective delay in the control loop may 
increase depending on the nature of the cause. This effective 
delay increase manifests itself as a phase shift between the 
commanded and actual system signals and can instigate PIOs. 
The proposed control allocator reduces the effective time delay 
by minimizing the phase shift between the commanded and the 
actual attitude accelerations. Simulation results are reported, 
which demonstrate phase shift minimization and recovery from 
PIOs. Conversion of the objective function to be minimized 
and constraints to a form that is suitable for implementation 
is given. 


I. Introduction 

A pilot induced oscillation (PIO) can be described as 
an inadvertent, sustained aircraft oscillation which is the 
consequence of an abnormal joint enterprise between the 
aircraft and the pilot [1]. There are several possible insti- 
gators of PIOs such as rate saturated actuators, high gain 
pilot/controller, system delays and phase lags. The focus 
of this paper is the rate saturation which is frequently 
observed during PIO events and has led to several crashes. 
An excellent overview of the effect of rate limiting in PIOs 
can be found in [2]. 

Figure 1 presents a basic model for a rate limited actuator 

[2] , where u is the input to the actuator and S is the actual 
actuator deflection. Without the rate limit, this dynamics is 
simply a first order lag, which is often used as an approx- 
imate model for actuators. Figure 2 shows time evolutions 
of input-output signals of such a rate saturated actuator [2], 

[3] , [4], where u c = u represents the pilot command. Gain 
reduction and an increase in effective time delay are two 
detrimental results of rate saturation, as seen from this figure. 



Fig. 1. Actuator model with rate saturation. 
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Fig. 2. Input u and output (5 of a rate saturated actuator. 


While gain reduction can not be avoided, there are several 
successful approaches in the literature that address elimi- 
nating the effective time delay: The main idea in [5], [6], 
[7], [8], [9] for eliminating the effective time delay is the 
differentiate-limit-integrate (DLI) approach. This simple yet 
effective method is implemented using a “software rate 
limiter” as shown in figure 3, where the software limiter 
is placed between the command signal and the input signal 
to the actuator. This method eliminates the effective time 
delay introduced by the rate saturation as seen in figure 4 
and hence the onset of a PIO can be avoided. 



Fig. 3. Software rate limiter. 



Fig. 4. Input u c and output 5 of a rate saturated actuator with a preceding 
software rate limiter. 


The DLI method proved very successful both in simulation 
and actual flight tests. It has, however, some deficiencies 
such as introduction of a bias and susceptibility to noise. 
See figure 5. These problems may be handled using different 
techniques such as filtering and resetting/retrimming [6], [7], 
[8]. See also [9] for more improved results. 



Fig. 5. Bias caused by the differentiate-limit-integrate approach 


There are also other methods for dealing with the effective 
time delay based on manipulating the input signal: In [10] 
a nonlinear adaptive filter is used to attenuate the pitch 
stick shaping function gain depending on the magnitude 
and frequency of the pilot input and the rate limits of the 
elevators. This approached was successfully implemented in 
Space Shuttle Enterprise ALT-5 and no further pitch PIO 
events were reported in the open literature since the 1977 
event [2]. In [11], a phase compensating filter is used which 
reverses the direction of a rate saturated actuator when the 
derivatives of the input and the output signals have opposite 
sign and when the absolute value of the error is increasing. 
In [12], a similar actuator output reversing logic is developed 
wherein a feedback signal is passed through a lowpass filter. 
This approach also proved to be successful in actual flight 
tests. In [13], another phase compensator is developed for 
a rate limiting element without using any feedback or logic 
but employing describing function relationships developed in 
[14]. This filter also performed successfully in flight tests. 

To the best of authors’ knowledge, all the previously 
reported successful implementation results were for SISO 
applications without any redundant actuators. Consider the 
closed loop flight control structure in figure 6, where the 
pilot is also in the loop. In this structure, the pilot task 
may be to track an altitude reference r, by getting altitude 
measurement feedback y, and making necessary corrections 
via a pilot stick which gives pitch rate commands u c as a 
reference to the inner flight controller. The inner controller 
may also be responding to roll and yaw rate commands 
at the same time. So, the “pilot command” u c can be a 
vector of three elements. The controller then calculates the 
necessary attitude accelerations v G 5ft 3 and then control 
allocator allocates the available actuators u E 3? m , m > 3, to 
achieve these desired accelerations while possibly satisfying 
secondary objectives like drag minimization. In this scenario, 
it is not obvious where and how to use the DLI software 
limiter. An extension of using the DLI method to multi-input 
multi-output (MIMO) applications is given in [6], however, 
the authors had to use “ganged” actuators for successful 
implementation. Ganging of the actuators, on the other 
hand, prevents the use of redundant actuators for secondary 
objectives like drag minimization or reconfiguration after a 
failure. In addition, ganging becomes more cumbersome as 
the number of actuators increases [15]. 

The novel control allocation method proposed in this paper 



Fig. 6. Overall SISO system structure 


builds upon the previous works and is suitable for MIMO ap- 
plications in the presence of redundant actuators. For ease of 
referencing, this technique will be called Control Allocation 
(to recover from) Pilot-Induced Oscillations (CAPIO). The 
main idea behind CAPIO, for which a patent is pending, 
is to minimize the phase lag introduced into the system 
due to rate saturation by minimizing the error between the 
derivatives of desired and actual total control effort vectors as 
well as minimizing the error between them, using constrained 
optimization techniques. To achieve this goal, for example in 
a SISO case, one needs to minimize the phase lag between 
the pilot input and the control surface deflection. On the 
other hand, in a MIMO case, where there are multiple inputs 
and outputs, one needs to pinpoint where exactly the phase 
lag is being introduced to the system. For example, in a 
scenario where the flight control system produces the desired 
rate accelerations and a control allocator distributes these 
commands to redundant actuators using some predefined 
optimization routine, it makes more sense to minimize the 
phase lag between the desired and achieved accelerations 
than concentrating on individual actuator signals. It is noted 
that merely having a control allocation scheme that takes 
into account the rate limits of the actuators as constraints 
can not prevent phase shift between the desired and achieved 
accelerations when saturation is unavoidable, and thus may 
not be able to handle a PIO situation. In simulation studies, 
where PIOs are present with conventional control allocation 
techniques and for high gain pilot models it will be shown 
that the onset of these PIOs can be prevented using CAPIO. 

It is noted that in [16] Durham and Bordignon extended 
the direct control allocation scheme to make it easier to 
implement for the case of rate-limited actuators and conse- 
quently ended up with a ’’moment-rate allocation” scheme. 
Although in [16] there is no implementation result showing a 
PIO preventing example, this control allocation scheme has 
a potential to handle PIOs despite being more complicated 
than CAPIO. Furthermore, the technique in [16] needs the 
calculation of a moment rate set which can introduce addi- 
tional computational intensity. 

The organization of this paper is as follows. In Section 
II, CAPIO is explained using a simple scalar example. In 
Section III, CAPIO is implemented for an aircraft model 
with redundant actuators, which is the case for which the 
use of CAPIO is particularly advantageous. In all cases, 
first, an example of a PIO formation is presented when a 
conventional control allocation is used, then it is shown that 
with the CAPIO the aircraft can recover from the PIO event. 
Finally, a summary of the paper is given. 





II. Scalar example 

In this section, it is assumed that the “controller” in 
figure 6 does not exist and hence, only the pilot is responsible 
for the control of the aircraft. The resulting system structure 
is presented in figure 7. In this scenario, the task is tracking 
a desired pitch angle Od . To simplify the analysis, the pilot 
is modeled as a static gain acting on the error between the 
desired pitch angle 6d and the measured pitch angle 6. It 
is noted that often pilot behavior is assumed to have at 
least a lead compensation while out of a PIO event, yet 
a pure gain during a fully developed PIO gives reasonable 
accuracy [3], [4], [17]. For our purposes, it does not make 
a difference whether a pure gain or a lead is used, except 
that the “gain only” option is simpler to work with. For 
an overview of different pilot models and behaviors see the 
lecture by McRuer [1]. The actuators used in this scenario are 
the elevators that are ganged. For the purpose of illustration, 
it is assumed that the elevators have an 0.05 second time 
constant with a rate limit of 28.7 deg/s and without any 
position limits. Short period approximation model [18] is 
used for the aircraft dynamics, which is given as 

0(s) _ 1.39(5 + 0.306) 

S e (s) ~ s 3 + 0.805s 2 + 1.3255 ( } 

where S e is the elevator deflection. It is noted that the elevator 
dynamics together with the rate limit is used in series with 
(1) in the simulations. 



Fig. 7. Overall SISO system structure 


A. Pitch angle control with conventional control allocation 

A control allocator distributes the total control effort, 
v E 5R P , to redundant actuators, u E 5ft m , m > p, so that 
total control effort is produced by the actuators through an 
input matrix, B E 5R pxm , possibly together with a secondary 
objective like drag minimization. One way of achieving 
this goal is using mixed optimization [15], which requires 
minimizing the following objective function 

J = || Bu - v\\ + e\\u - u p \\t (2) 

subject to U'TTT.m u ^ u max ^nd Umin ^ u ^ Umax, where 
u p is a “preferred” actuator position and e is a weight. This 
weight may be chosen sufficiently small to emphasize the 
main goal of minimizing the error between the desired total 
control effort v and the produced control effort Bu by the 
actuators. 

In our example, the total control effort, u c = v, is a scalar 
and there is only one effective actuator (ganged elevators), 
which makes B = 1. To simplify further, it is assumed 
that there are no secondary objectives, thus 6 = 0 and 


the elevators have only rate limits. As a result, the control 
allocator needs to minimize 


J = \u — u c \, (3) 

subject to iimin < u < Umax • In this form, the main task of 
the control allocator becomes trying to achieve u c as close as 
possible while preventing the rate saturation. However, it is 
noted that the control allocator will still make the elevators 
work as close as possible to rate saturation limit to be able to 
follow u c , whenever necessary. Therefore, one will observe 
almost the same input output relationship as in figure 2 
between u c and S e . As a result, it turns out that a conventional 
control allocator for a scalar system is equivalent to a pure 
gain of 1, in terms of the system response. 

Figure 8 presents a simulation result where a 1 radian 
pitch step is given as a reference to the pilot and the pilot 
acts on the tracking error with a gain of 1.2. Although during 
the first 4.5 seconds the actuator is rate saturated, after 4.5 
seconds the pilot command (u c ) rate gets smaller and thus 
the actual elevator deflection S e can catch it up. The result 
is a stable tracking of the step reference. 



Fig. 8. Pitch reference tracking with conventional control allocation 


Figure 9 presents a simulation result when a more aggres- 
sive pilot with a gain of 1.65 is controlling the aircraft. Due 
to a more aggressive pilot control, the elevators are always 
saturated and an introduced effective time delay is visible 
(note the similarity with figure 2). The result is formation 
of a PIO observed as sustained oscillations in the system 
response. 

Figure 10 presents the tracking error and the pitch angle 
in the same plot. It is noted that the signals have the same 
magnitude with a —180 degrees phase shift, which shows 
that the closed loop system is pushed to the marginally stable 
point with the introduction of an effective time delay. This, 
in fact, confirms that this is a PIO event, not just temporary 
maneuvers conducted by the pilot to stabilize the system. As 
in a typical PIO event, the pilot is not driving the oscillation 
but she is driven by the oscillation, meaning that she must 
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Fig. 9. PIO formation during pitch reference tracking with conventional 
control allocation. 
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Fig. 11. Bode plot of the overall system without the rate limiting element. 


redirect efforts away from the primary task by a noticable 
amount [17]. 



Fig. 10. Evolution of pitch angle and pitch angle tracking error. 


Figure 1 1 shows the bode plot of the loop transfer function 
neglecting the rate limiting element. The phase margin is 
23.2 degrees with a crossover frequency of 1.77 rad/s which 
gives a 0.228 seconds time delay margin. Introduced effec- 
tive time delay is larger than this value (see figure 9) which 
actually makes the system unstable (instead of marginally 
stable). It is also noted that there is a gain reduction due to 
rate limiting which increases the time-delay margin. 

B. Pitch angle control with CAPIO 

As discussed in the previous section, the effective time 
delay introduced due to rate saturation can be given as one 
of the reasons for the formation of the PIO observed in 
figure 9. This effective time-delay is the result of the phase 
lag between u c and S e , where “phase lag” refers to the 
shift between the extremum points of the signals. Hence, the 
main goal of CAPIO is to minimize the phase lag between 


these two signals. This is achieved by modifying the original 
objective function (3) as 


J = \u — u c \ + k \u — u c | , (4) 

where k > 0. It is noted that with this modified objective 
function, the control allocator is trying to realize ii c as well 
as u c . A very high value of k makes the two signals, u c 
and S e , have approximately the same derivative at all times, 
which eliminates the phase lag completely but causes a bias 
as shown in figure 5. On the other hand, a very small value 
of k may not be sufficient for the control allocator to be 
any different than the conventional one and thus does not 
prevent PIOs. Therefore, the designer needs to decide on a 
suitable value of k which minimizes the phase lag and at the 
same time prevents a bias. As an alternative, the designer 
can choose to “activate” k, i.e. set it to a positive value, 
only when it is needed, and keep it 0 at all other times. The 
latter approach is taken in this paper, assuming that there 
exists a PIO detection algorithm on board. 

For implementation, the objective function (4) needs to 
be put in a form that can be minimized using available 
optimization packages. To achieve this goal, the derivatives 
in the objective function are approximated as follows 


u(t) » “W-"(>- r ) 


( 5 ) 


where T is the sampling interval. In all the simulations in 
this paper, we use T = 10ms. Using the approximation in 
(5), the following objective function is obtained 


J (t) = T\u(t) — u c (t)\ + k\u(t) — u(t — T) — Tu c (t ) | (6) 


where, j' = TJ. 

Figure 12 shows the simulation result where the system 
is provided with the same pitch angle reference and the 
same pilot gain as in figure 9, but this time CAPIO is the 
governing control allocator, k > 0 is set with the detection 
of the PIO and deactivated (set to 0) when the system is 


out of the PIO. With the activation of k , actual elevator 
deflection 5 e comes in phase with the commanded elevator 
deflection u c and this eliminates the previously introduced 
effective time-delay and drives the system out of the PIO. 
Then, k is deactivated and CAPIO continues to perform as 
a conventional control allocator ensuring that 5 e tracks u c . 
As a result, the system eliminates the PIO without any bias 
formation. 



Fig. 12. Pitch reference tracking with CAPIO 


Remark 1 The objective function (6) contains the 
derivative of the commanded elevator deflection u c . This 
derivative can be computed using differentiation filters, such 
as s/(ts + 1). The measurement noise can be alleviated 
by using low pass filters, which proved to be successful in 
previous PIO preventing algorithms [9]. 


are redundant actuators in the system. 


III. Multi Input Multi Output (MIMO) Example 
with Redundant Actuators 

To illustrate the advantages of CAPIO, a flight control 
example using a simplified ADMIRE model from [19] is 
used. This model includes redundant actuators which makes 
the DLI method hard to apply if one does not want to gang 
the actuators. 

The linearized aircraft model at Mach 0.22, altitude 3000m 
is given by 


x = [a /? p q r] T -x i in , 
y = Cx=\p q r] T — y Un , 

S = [<5 C S re Sle Sr] 1 ' — (5ii„, 

U = [u c U re Ule U r ] T - Ulin 
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where a, /3, p, q and r are the angle of attack, sideslip 
angle, roll rate, pitch rate and yaw rate, respectively. S and 
u represent the actual and the commanded control surface 
deflections, respectively. Control surfaces are canard wings, 
right and left elevons and the rudder. (.)i in refers to values at 
the operating points where the linearization was performed. 
The actuators have the following position and rate limits 

Sc e [-55, 25] x 22; 6 re , S u , S r € [-30, 30] x 22 

S c ,6 re ,6 le ,6 r e[- 70, 70] x -2- (8) 


Remark 2: PIO detection algorithm that is assumed 
to be present on board can detect the PIO earlier or later 
than what is assumed in figure 12. In addition, this algorithm 
can decide that the system is out of the PIO earlier or later 
than what is assumed in figure 9. These timing decisions 
can affect system performance. However, the main point 
made here is that although the timing of the activation and 
deactivation of the k term may effect the overall system 
performance, once k is activated, CAPIO will force the 
actuator deflection to be in phase with the commanded 
signal and eliminate the introduced effective time delay. 
The CAPIO will continue to perform as a regular control 
allocator once k is deactivated. This dual behavior both 
helps the aircraft go out of the PIO and prevents bias 
formation. 

Remark 3 This scalar example is given mainly for a 
simple introduction to the CAPIO. It is clear that the 
DLI method would also prevent the PIO in this scenario. 
However, the real power of the CAPIO becomes apparent in 
more realistic cases where the system is MIMO and there 


and have first order dynamics with a time constant of 0.05 
seconds. It is noted that the position limits given are the same 
as in [19] but the rate limits are assumed to illustrate CAPIO 
properties. At trimmed flight, we have 

x = x 0 = (12.7 0 0 0 0) 

5 = S 0 = (0 5.4 5.4 0) (9) 

To make this model suitable for control allocation im- 
plementation, the actuator dynamics are neglected and the 
control surfaces are viewed as pure moment generators and 
their influence on a and /? is neglected. It is noted that the 
actuators dynamics are present during the simulations, i.e. 
they are neglected only during the control allocation algo- 
rithm derivation. These assumptions lead to the following 
approximate model 

x = Ax + B u u = Ax + B v v , 
v = Bu , 


( 10 ) 


where 


B u = B V B , 
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- 0.1179 
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0 

- 0.9661 

0 

- 10.5128 

- 0.9967 

0 
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2.6221 
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0 

- 0.5057 
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B = 


0 0.7075 - 0.0939 0 - 0.2127 

0 -4.2423 4.2423 1.4871 

1.6532 -1.2735 -1.2735 0.0024 

0 -0.2805 0.2805 -0.8823 


The virtual (total) control effort, v, consists of the angular 
accelerations in roll, pitch and yaw. 

In this flight control example the pilot task is to track a 
given pitch angle reference, Qd , using a pitch rate, qd, stick. In 
addition, roll rate, p , and the yaw rate, r, are to be controlled 
independently to track their references pd and r^. The overall 
system structure is given in figure 13. 
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Fig. 13. Overall MIMO system structure 


A. Flight control with conventional control allocation 

The conventional control allocation used in this example 
minimizes the following objective function 

J=\\Bu-v\\ 2 2 + e\\u\\l (13) 

subject to max(u min T + u k -i,u min ) <u< mm(u max T + 
Uk-hUmax), where T is the sampling interval. It is noted 
that norms, instead of square-norms, can be used in the 
objective function. Note that (13) is in the form of a typical 
objective function used in conventional control allocators 
[15], where the main objective is to minimize the error 
between the desired and the actual total control efforts. As 
e — > 0, minimizing (13) becomes equivalent to achieving 
the main objective explained above and picking the solution 
that gives the minimum control surface deflection, among 
different solutions. In this example e = 10 -5 . 

Figure 14 presents the simulation result with the conven- 
tional control allocation where the pilot receives a step pitch 
angle reference at t = 3 seconds and the inner controller 
receives a pulse roll rate reference at t = 0.5 seconds 
and a zero yaw rate reference at all times. The pilot is 
aggressive and has a gain of 4.07. Because of this high gain, 
the aircraft goes into a divergent PIO in the pitch axis. In 
addition, yaw axis also goes into a sustained oscillation mode 
which diverges towards the end. The inner controller can still 
manages to track the roll rate pulse reference. Canard wings 
and the ailerons saturate both in position and the rate. The 
results of this saturation can best be observed as a phase shift 
between the desired pitch acceleration v 2 and the actual pitch 
acceleration created by the control surfaces Bu<i , although 
the same phenomenon is observed between i/ 3 and Bus. 
These phase shifts, or the effective time delays, is something 
that is almost always observed in PIO events due to actuator 
saturation. 


The inner controller is a dynamic inversion controller 
which uses qd , Pd and as references and produces the 
necessary attitude accelerations, v G 5P 3 , to track these 
references. Dynamic inversion control laws, v, make the 
closed loop dynamics follow a desired reference model 

Vm = AmUm H - RmXm Of) 

where y m = [p m q m r m ] T represents the desired output 
and r m = [pd qd r d] T is the reference input. In this 
example, A m = —2 x / 3x3 and B m = 2 x J 3x 3. Reference 
model tracking can be achieved by inverting the dynamics 
[15] as 

v = (CB v )- 1 [A m y + B m r m -CAx}. (12) 

The control allocator distributes this total control effort, 
v, to individual control surfaces via the actuator commands, 
u G 5ft 4 . The control surfaces then produces actual attitude 
accelerations, Bu, where B is the control input matrix. The 
pilot is modeled as a pure gain for simplicity. 


B. Flight control with CAPIO 

To recover from a PIO event, CAPIO forces the virtual 
(total) control effort v, to be in phase with the actual control 
effort Bu produced by the actuators. To achieve this, a 
derivative error term is added to objective function (13) to 
obtain the following CAPIO objective function 

/ = || Bu- v\\ 2 2 + || W d (Bu - v)\\\ + e||u||l (14) 

where Wd G R 3xS represents a weighting matrix on the 
derivative term. The cost function J is minimized with 
respect to u, with u = (u — u~)/T, where u~ denotes the 
value of u at the previous sampling instant. As discussed 
in the previous section, one can affect the phase lag by 
changing the weighing matrix. In addition, one can achieve 
axis prioritization by properly choosing Wd . As in the scalar 
case, it is assumed that there exits a PIO detector on board 
and we use an implementation in which Wd is set to identity 
or zero depending on whether the aircraft is in a PIO event 
or not. 

The objective function (14) needs to be transformed into a 
form that can be minimized numerically. To achieve this goal, 







Fig. 14. Pitch angle 0 and aircraft states x , on the left. Desired 
(commanded) and actual attitude accelerations v and Bu, and the control 
surface deflections 5, on the right, when a conventional control allocator is 
used. 


Fig. 15. Pitch angle 0 and aircraft states x , on the left. Desired 
(commanded) and actual attitude accelerations v and Bu, and the control 
surface deflections (5, on the right, when CAPIO is used. 


the derivatives in the objective function are approximated as 
in (5). After some algebra (14) can be rewritten as 

j = u T ( B t T 2 B + B t RB + 6 / 4 x 4 ) U 

+2 (- v t T 2 B - u~ T B T RB - v t TRB ) u 

+v t T 2 v + u~ T B T RBu~ + 2 u~ T B t RTv 
+v t T 2 Rv (15) 

subject to max(u min T + u~,u min ) < u < min (u max T + 
u~,Umax), where R = W d T W d . 

Figure 15 presents the simulation result when CAPIO is 
used as the control allocator. All the settings including the 
pilot gain are the same as in the previous example with the 
conventional control allocation. Since CAPIO prevents the 
effective time delay introduction, the aircraft now recovers 
from the PIO. 

To show the difference that CAPIO makes in control 
effort realization, the pitch axis accelerations are presented 
again in figure 16 for both cases. It is noted as as soon 
as the PIO detection signal is obtained, CAPIO forces the 
control surfaces to produce accelerations “in phase” with 
the commanded accelerations, eliminating the time delay 
due to phase shift. When the aircraft recovers from PIO, 
control surfaces arrange themselves to track the commanded 
acceleration. The result is recovery from the PIO without any 
bias formation. 


The constrained optimization of the cost (15) is a low 
dimensional quadratic programming problem with linear 
inequality constraints. This problem depends on a vector 
of parameters, specifically, on p = [u~ v v\. Note that 
the parameters enter linear in the cost and in the constraints 
and, hence, such a quadratic programming problem can be 
solved explicitly using off-line multi-parametric QP solvers. 
The solution is known to be a piecewise affine continuous 
function of the parameter vector and have the following form, 
u = fip + gi, if Fi + Gi < 0, i = 1, ..., N r where N r is 
the finite number of polyhedral regions and each region is 
associated with its set of linear inequality constraints and its 
affine map. Such an explicit solution is computed off-line 
and can be deployed on-line in the aircraft software using a 
set of simple if-then-else rules, additions, multiplications and 
comparisons. The need to embed a quadratic programming 
solver to perform constrained optimization of the cost (15) 
within aircraft software can thus be avoided altogether. A 
cross-section of the explicit solution polyhedral regions is 
illustrated in figure 17; the explicit solution has N r = 223 
regions. 

IV. Summary 

In this paper, a new control allocation scheme for recov- 
ery from pilot induced oscillations (CAPIO) was proposed. 
CAPIO functions by minimizing the error between the 
derivatives of desired total control effort and the achieved 



Fig. 16. Pitch angle 0 and aircraft states x , on the left, desired 
(commanded) and actual attitude accelerations v and Bu, and the control 
surfaces S, on the right, when a CAPIO is used. 


Polyhedral partition - 223 regions 



v 2 =pitch acceleration 

Fig. 17. Cross-section of regions of explicit solution by pitch acceleration 
- pitch acceleration rate plane. 


total control effort vectors, as well as minimizing the error 
between these two vectors. The derivative error minimization 
forces the elements of these two vectors to have the same 
phase and the error minimization makes them converge to 
each other. This dual behavior results in recovery from 
dangerous PIO events without creating any bias. The tech- 
nique was explained first on a scalar example, and then 
a MIMO example with redundant actuators was given. In 
each example, first a PIO was created using high gain 
pilot models with a conventional control allocator and then 
recovery from this PIO event was shown by the employment 
of CAPIO. Conversion of the allocation scheme to a form 
that is easy to implement was also given. It is noted that 
in real aircraft implementation, handling of measurement 


noise, digital realization of the derivative and the integration 
of CAPIO with a PIO detector are important points to be 
addressed for a successful technology transfer. 
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